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Abstract 

We have calculated the next-to-leading order cross sections for the inclusive produc- 
tion of L'*-mesons in ep collisions at HERA for finite, although very small Q^. In 
this Q^-range, the same approximations as for photoproduction can be used. Our 
calculation is performed in the general-mass variable-flavour-number scheme. In 
this approach, large logarithms of the charm transverse momentum are resummed 
and finite terms depending on m? /pj, are kept in the hard scattering cross sections. 
The theoretical results are compared with recent data from the ZEUS collaboration 
at HERA. On average, we find good agreement. 
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1 Introduction 



The understanding of the dynamics of charm quark production at HERA has been im- 
proved considerably over the last ten years by the HI and ZEUS collaborations, who have 
performed many measurements of inclusive D*-meson production in the photoproduction 
mode of ep collisions with almost vanishing virtuality (Q^ ^ 0) of the exchanged pho- 
ton, as well as in the deep- inelastic (DIS) mode (with photon virtuality > 0). The 
theoretical description of heavy quark production in the framework of perturbative QCD 
is complicated due to the presence of several large scales appearing in this process. In 
DIS, is large, but also the transverse momentum pj- of the produced D*-meson may 
be large. In addition, depending on the kinematic range considered, also the mass m 
of the charm quark may have to be taken into account. Different calculational schemes 
have been developed which can be applied for an interpretation of experimental data, 
depending on the specific kinematical region and the relative importance of these three 
scales. 

In the case of relatively small transverse momentum, px ^ m, the fixed-flavour number 
scheme (FFNS) is usually applied. Here one assumes that the light quarks {u, d, s) and 
the gluon are the only active flavours within the proton and the photon. In this scheme 
cross sections for e + p ^ e' + D* + X have been calculated for DIS in Ref. [1] and for 
photoproduction in Ref. [2]. In photoproduction, where ~ 0, the direct process has to 
be supplemented with the resolved process, where the photon acts as a source of partons 
which interact with partons in the proton. These two interaction modes are needed to 
describe the singular region at = for massless quarks appropriately. In the FFN 
scheme [I] E] the charm quark appears only in the final state of the direct and resolved 
processes, via the hard scattering of light partons, including the photon. The charm quark 
mass m is explicitly taken into account together with the transverse momentum of the 
produced D*-meson; this approach is therefore expected to be reliable when p^ and m 
are of the same order of magnitude. In the FFNS, the charm quark mass acts as a cutoff 
for the initial- and final-state coUinear singularities and sets the scale for the perturbative 
calculation. The mass m is fully retained in the calculation of the hard-scattering cross 
sections. 

In the complementary kinematical region where pt ^ m, calculations are usually based 
on the zero-mass variable- flavour- number scheme (ZM-VFNS). This is the conventional 
parton model approach where the zero-mass parton approximation is applied also to the 
charm quark, although its mass is not small and large compared with Aqcd- In the ZM- 
VFNS, the charm quark acts also as an incoming parton with its own parton distribution 
function (PDF) in the proton and in the photon leading to additional direct and resolved 
contributions. Usually, charm quark PDFs and also the fragmentation functions (FFs), 
describing the transition c D*, are defined at an initial scale /ig chosen equal to the 
charm mass m. Then this is the only place, where the charm mass enters in this scheme. 
The D*-meson is produced by fragmentation from the charm quark produced in the hard 
scattering process; but also fragmentation from the light quarks and the gluon has to be 
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taken into account. The well-known factorization theorem provides a unique procedure 
for incorporating the FFs into the lowest order (LO) and next-to-leading-order (NLO) 
perturbative calculations. The predictions obtained in this scheme are expected to be 
reliable only in the region of large pt since all terms of the order rn^ /p\ are neglected in 
the hard scattering cross section. Calculations for D*-production in the ZM-VFNS have 
been performed some time ago for photoproduction in Ref. [3j and for DIS in Ref. 

A unified scheme that combines the virtues of the FFNS and the ZM-VFNS is the so- 
called general-mass variable-flavour-number scheme (GM-VFNS) [3]. In this approach 
the large logarithms ln(p|,/m^), which appear due to the coUinear mass singularities in 
the initial and final state, are factorized into the PDFs and FFs and summed by the 
well known DGLAP evolution equations [6J for the PDFs and FFs. The factorization is 
performed following the usual MS prescription which guarantees the universality of both 
PDFs and FFs. At the same time, mass-dependent power corrections are retained in the 
hard-scattering cross sections, as in the FFNS. It is expected that this scheme is valid not 
only in the region ^ m?, but also in the kinematic region where pt is larger than a 
few times the charm mass m only. We should emphasize that in the GM-VFN scheme, 
the incorporation of the fragmentation c D* is based on the factorization theorem; 
this is a prerequisite for applying the FFs in different processes. In the usual FFNS 
calculation, this is not the case and non-perturbative FFs for the transition c ^ D* can 
be supplemented on purely phenomenological grounds only. 

It is the purpose of this work to present theoretical results for the D*-production cross 
section in e^p scattering with finite non-zero photon virtuality in the region 0.05 < < 
0.7 GeV^ and discuss a comparison with experimental results obtained with the ZEUS 
detector at HERA We shall calculate all cross sections with the same kinematical 
constraints as in the ZEUS analysis [7j in the GM-VFNS. Since in the ZEUS experiment 
the photon virtuality is small, the application of the photoproduction approximation is 
justified where the Q^-dependence in the hard scattering cross sections is neglected. 

Details of our calculation have been described recently in Ref. |H]. In this work, we had also 
studied theoretical uncertainties of the photoproduction cross section in the GM-VFNS 
due to various possible choices for input variables, as for example, the proton and photon 
PDFs and the D* FFs. In that reference one can also find a discussion of the dependence 
on the factorization and renormalization scales and the influence of the charm quark mass 
m. The application of the theoretical framework described there to the present case of 
low-Q^ DIS is straightforward and amounts to an adjustment of the parameters entering 
the Weizsacker- Williams approximation for the flux of the virtual photon. Instead of 
Qmin = ^ly'^/i^ ~ y) cLnd the value for Q^ax given by the anti-tagging condition of the 
final electron (positron) as used in the measurement of the photoproduction process, we 
have to fix Q^m Q^ax ^^e values used in the ZEUS analysis. We will consider the 
kinematic range as in Ref. 0, i.e. 0.05 < Q'^ < 0.7 GeV^. 

The outline of the paper is as follows. In Section 2 we give a short description of the 
various input options used in the calculation. Section 3 contains our results and the 
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comparison with the experimental data from ZEUS. 



2 Input choices for the calculation 

The D*-electroproduction cross section o"ep(v^) at the ep center-of-mass energy ^/s is 
related to the photoproduction cross section at center-of-mass energy W^p, cr-ypiW^p), in 
the following way: 

{yV~s) . (1) 



mm 



Here, f^-y is the energy spectrum of the exchanged virtual photon which in the Weizsacker- 
Williams approximation [9j is given by 



y ^min \ 'Vm 




(2) 



The photon flux /e^ depends on y, Ql^^^ and Q^^^. The range of y, y^m < y < 1/max, 
as well as the limits Qmin Qmax^ ^'^^ determined by the cuts and bin limits in the 
experimental analysis. In photoproduction, Q^in ^ is very small, whereas in the 
ZEUS analysis Qmin = 0-05 GeV^ and Q^ax = 0-''' GeV^, or within these limits for the 
measurement of the distribution, a is the electromagnetic fine structure constant and 
y = E^/Eg, the ratio of the energies of the incoming photon and electron, is determined 
by the inelasticity y = Q"^/ {2P ■ q) where P and q are the 4-momenta of the incoming 
proton and the photon. 

The cross section for direct photoproduction in Eq. ([1]) is a convolution of the proton 
PDF, the fragmentation function for the transition of parton a to the observed Z}*-meson 
(where a = u, u, d, d, s, s, c, c, and g) and the cross section for the hard scattering process 
76 —>■ aX. For the resolved contribution, an additional convolution with the photon PDFs 
has to be performed. The hard scattering cross sections are calculated including next- 
to-leading order corrections of the order 0(a<j). The PDFs are evolved at NLO. For the 
photon PDF we use the set GRV92 of Ref. ^Q\, converted to the MS factorization scheme; 
for the proton PDF we have chosen the most recent parametrization CTEQ6.6M [llj of 
the CTEQ group. 

For the FFs we use the set Belle/CLEO-GM of Ref. [I2]. Note that in our previous 
work [S] we had used the set Global-GM instead, where also LEPl data [13] at large s 
had been included in the fit, whereas the set Belle/CLEO-GM is based on a fit of the 
combined Belle [H] and CLEO [15] data at -y/i = 10.52 GeV only. For the photoproduc- 
tion cross section da/dpT we found results larger by 25 — 30% in average when using the 
Belle/CLEO-GM parametrization, as compared to the set Global-GM of Ref. [12]. The 
strong coupling constant a^s^^\nji) is evaluated with the two- loop formula [TB] with rij = 4 
active quark flavours and the asymptotic scale parameter = 328 MeV, corresponding 
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to af'\mz) = 0.118. The charm quark mass is fixed to m = 1.5 GeV. We choose the 
renormahzation scale fiR and the factorization scales fip and fip/ related to initial- and 
final-state singularities to be fiR = (.rttit and = hf' = ^pf^T, where mT = m? + pj, 
is the transverse mass and and are parameters varied about their default values 
^j? = 'C-F = 1 in order to assess theoretical uncertainties as described below. 

As already mentioned in the introduction, the photoproduction cross section is calculated 
in the GM-VFN scheme. In this scheme the cross section has a smooth limit for m — >■ 
and approaches the result of the ZM-VFN scheme for m or — > oo. The basic 
features of the GM-VFNS are described in Ref. [8] and the literature quoted therein. 
Compared to the ZM-VFNS, the GM-VFNS incorporates mass-dependent terms as in the 
FFNS, so that the cross sections calculated in the GM-VFNS are supposed to be valid also 
for medium scale pt values, close to the heavy-quark mass. In order to conform with the 
MS factorization of singularities, finite subtraction terms must be supplemented to the 
results of the FFNS calculation. These subtraction terms had been calculated in Ref. [I7j 
for the direct photon, and in Ref. [TH] for the resolved photon contributions. They include 
logarithmic, scale-dependent contributions related to gluon emission from charm quarks 
and to charm- ant icharm production from incoming gluons. As a consequence, in the GM- 
VFNS one has to take into account also processes with incoming charm quarks, involving 
corresponding charm-quark components in the PDFs of the photon and the proton. In 
addition, for the final state, apart from the FF describing the transition c ^ D*, also 
FFs for the transition of a light parton to the heavy meson, a — > D*, are needed. These 
contributions are not present in the FFNS calculation. Instead, they are taken into 
account at fixed order of perturbation theory as part of the hard scattering cross sections. 
Compared with the ZM- and GM-VFNS approaches, the FFNS has the advantage that 
it is valid also for < px ^ m, a. property which is not realized in the presently available 
implementation of the GM-VFNS. A reliable prediction down to px = is necessary if 
total cross sections for heavy-meson production are to be calculated. 

In our calculation we implement the experimental conditions of the ZEUS analysis: the 
energies of the incoming protons and electrons (positrons) are Ep = 920 GeV and Eg = 
27.5 GeV, respectively. The inelasticity y varies in the range 0.02 < y < 0.85. The 
transverse momentum pt and the rapidity f] of the D*-meson have been measured in the 
kinematic ranges 1.5 < px < 9.0 GeV and \ri\ < 1.5. The photon virtuality is taken in the 
interval 0.05 < < 0.7 GeV^. 



3 Results 

In this section we present our results for the differential cross sections of the process e+p — » 
e'+D*+X in the GM-VFNS as a function of Q^, y, px and rj. We choose bins in accordance 
with the ZEUS measurement and differential cross sections are obtained by averaging over 
the bin sizes as are the experimental ones [7] . The numerical results are shown in Fig. [H 
We have estimated theoretical uncertainties by varying independently the parameters C,r 
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Figure 1: Differential cross sections for D*-meson production in low-Q^ ep scattering, 
compared with experimental results from the ZEUS collaboration [7]. The default choice 
for the scale parameters is C,r = C,f = ^ (full lines) and error bands are obtained by varying 
and (see text, dashed lines). The kinematic range is given by 0.05 < < 0.7 GeV^, 
0.02 <y< 0.85, 1.5 < pr < 9.0 GeV and \r]\ < 1.5. 



and in the range 0.5 < ^r, ■Cf < 2 about their default values C,r = ^f = ^ imposing the 
constraint 0.5 < ^f/^_r < 2. The maximal and minimal differential cross sections obtained 
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this way are shown in Fig. [T]as dashed hnes together with the central values (full lines). 
Comparing the errors of the ZEUS data points and the theoretical errors due to the 
scale variation we observe that all data points for da/dQ"^, da/dy, da/dpx and da/drj are 
compatible with our theoretical predictions. Most of the predictions with the default scale 
choice C,R = = ^ are found inside the range given by the data and their experimental 
errors. The discrepancies between theoretical prediction and experimental measurement 
are largest for the smallest values of Q^, px and y\ however, also the theoretical error 
band is largest for these bins. Since small values of these kinematic variables dominate 
at all values of 77, the uncertainty due to scale variations is large for all bins of da/drj. 
For da/dpT the relative size of the error band decreases with increasing px, as expected. 
In the first bin of the pt distribution where pT{mm) ~ 1.5 GeV, the scale variable fip 
is small, approximately equal to 1.5 x m, which explains the large sensitivity to scale 
variations. 

Actually we can not expect our theoretical framework to be particularly accurate at small 
Pt, close to px — fn, since the cross section in the GM-VFNS contains parts which are 
calculated in the massless approximation. This prevents us to calculate cross sections 
down to Pt = and consequently also the total D*-production cross section integrated 
over Pt can not be calculated reliably. In the kinematic region pt ~ m, the FFNS is 
the better choice. The comparison of the data in the first pt bin with the FMNR result 
for photoproduction [2], as presented in Ref. [7], shows indeed good agreement. Since 
da/dpT is over-estimated in our approach for pt — m, also the cross section integrated 
over the full range of pt values considered in the ZEUS experimental analysis is too 
large: the measured total cross section quoted in Ref. [7] i£] a{ep — ^ e'D*X) = lO.ltJJ 
nb, whereas the corresponding theoretical prediction is 13.9l']5°g'' ^b. In the restricted Ar- 
range 2.5 < < 9 GeV, we find instead 4.llo!9 nb to be compared with the experimental 
value of 3.9 ±0.4 nb (obtained from the results given in [7] by summing over the last three 
Pt bins), i.e. for larger pt the theoretical cross section in the GM-VFNS agrees with the 
measured one quite well. 

In obtaining theoretical predictions for the photoproduction cross section, one has to 
distinguish direct and resolved components. Each of the two parts depends on the choice 
of a factorization scheme and can not be compared separately to experimental data. 
Only the sum of the direct and resolved contributions has a physical meaning. From the 
theoretical point of view, it is nevertheless interesting to study the decomposition of the 
cross section into these two components, in particular in view of a comparison with results 
obtained in the FFNS approach for photoproduction [2] or for DIS pj. In DIS at large 
Q^, there is only a direct contribution to the cross section since initial-state singularities 
appearing in the limit ^ have not to be subtracted. The photoproduction cross 
section in the FFN scheme [2] has a smaller resolved part than in the GM-VFNS, since 
it originates only from initial gluons and light quarks, whereas in the GM-VFNS also 
processes with charm quarks in the initial state contribute. Therefore, in our case the 

^Statistical and systematic errors, and an error due to the uncertainty in the branching ratios, are 
added in quadrature. 
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Figure 2: Differential cross sections for D*-meson production in low-Q^ ep scattering in 
tlie kinematic range 0.05 < < o.7 GeV^, 0.02 < y < 0.85, 1.5 < pr < 9.0 GeV, 
|?7| < 1.5. The full cross section (full lines) is separated into its direct (dotted lines) and 
resolved parts (dashed lines) for C-R = = 1- 



direct and resolved contributions are of comparable magnitude. This is seen in detail in 
Fig. [2] where we have plotted these two components and their sum for da/dQ"^, da/dy, 
da/dpT and da/drj. 
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We conclude that in the general-mass variable-flavour-number scheme, the photoproduc- 
tion approximation, i.e. the approximation where the hard scattering cross sections are 
evaluated with = 0, leads to theoretical predictions in very good agreement with re- 
cent ZEUS data for finite, but small 7^ 0, in particular for the larger values of pt- The 
description of the measurement at the smallest values of px — 1-5 GeV is less satisfactory, 
but still in agreement with data within errors. 
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